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Uniqueness  in  remote  sensing  of  the  inherent 
optical  properties  of  ocean  water 


Michael  Sydor,  Richard  W.  Gouid,  Robert  A.  Arnone,  Vladimir  I.  Haltrin,  and  Wesley  Goode 


We  examine  the  problem  of  uniqueness  in  the  relationship  between  the  remote-sensing  reflectance  (Rrs) 
and  the  inherent  optical  properties  (lOPs)  of  ocean  water.  The  results  point  to  the  fact  that  diffuse 
reflectance  of  plane  irradiance  from  ocean  water  is  inherently  ambiguous.  Furthermore,  in  the  400  < 
X  <  750  nm  region  of  the  spectrum,  Rrs(X)  also  suffers  from  ambiguity  caused  by  the  similarity  in 
wavelength  dependence  of  the  coefficients  of  absorption  by  particulate  matter  and  of  absorption  by  colored 
dissolved  organic  matter.  The  absorption  coefficients  have  overlapping  exponential  responses,  which 
lead  to  the  fact  that  more  than  one  combination  of  lOPs  can  produce  nearly  the  same  Rrs  spectrum.  This 
ambiguity  in  absorption  parameters  demands  that  we  identify  the  regions  of  the  Rrs  spectrum  where  we 
can  isolate  the  effects  that  are  due  only  to  scattering  by  particulates  and  to  absorption  by  pure  water. 
The  results  indicate  that  the  spectral  shape  of  the  absorption  coefficient  of  phytoplankton,  aph(X),  cannot 
be  derived  from  a  multiparameter  fit  to  Rrs(X).  However,  the  magnitude  and  the  spectral  dependence 
of  the  absorption  coefficient  can  be  estimated  from  the  difference  between  the  measured  Rrs(X)  and  the 
best  fit  to  Rrs(X)  in  terms  of  lOPs  that  exclude  aph(X).  ©  2004  Optical  Society  of  America 
OCIS  codes:  010.4450,  010.7340,  290.4210,  290.7050. 


1.  Introduction 

This  investigation  was  prompted  by  two  questions 
raised  at  the  recent  Ocean  Optics  Conference  in 
Santa  Fe,  New  Mexico,  on  18~23  November  2002: 

1.  Why  do  we  have  trouble  with  remote-sensing 
algorithms  for  the  inherent  optical  properties  (lOPs) 
of  ocean  water  in  certain  areas  of  the  oceans  where 
colored  dissolved  organic  matter  (CDOM)  is  present; 
i.e.,  why  are  algorithms  for  lOPs  not  global? 

2.  \^y  do  we  need  hyperspectral  measurements 
of  remote  sensing  reflectance  [Rrs(X)]  rather  than 
measure  Rrs  at  few  well-placed  wavelengths  to  re¬ 
trieve  lOPs  by  use  of  remote-sensing  data? 

These  questions  lead  us  to  examine  the  uniqueness  of 
the  spectral  shapes  of  Rrs(X)  as  a  fimction  of  wave¬ 
length  X.  Mobley^  raised  the  question  of  uniqueness 
some  decade  ago.  Uniqueness  implies  that  we  can- 
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not  have  two  different  sets  of  parameters  that  predict 
the  same  reflectance  spectrum  within  the  experimen¬ 
tal  accuracy.  To  test  uniqueness  or  to  prove  the 
ambiguity  of  Rrs(X)  we  need  to  examine  its  mathe¬ 
matical  definition  and  devise  two  or  more  distinct 
sets  of  parameters  or  two  mathematical  solutions 
that  5deld  equivalent  Rrs  spectra. 

To  avoid  possible  preselection  of  data  in  testing  for 
ambiguity  of  Rrs(X)  we  examine  data  published  by 
other  investigators  as  well  as  our  own  recent  data 
from  experiments  from  the  Gulf  of  Mexico.  In  par¬ 
ticular,  we  make  impromptu  use  of  data  presented  by 
Roesler  and  Boss^  at  the  scientific  conference  men¬ 
tioned  above.  The  examination  of  Rrs  from  widely 
disparate  sources  and  geographic  regions  will  lend 
global  character  to  our  investigation  and  will  leads  us 
to  a  better  imderstanding  of  the  underlying  physical 
processes  that  give  rise  to  Rrs. 

Let  us  first  examine  the  accepted  mathematical 
expressions  for  diffuse  reflectance  R  from  ocean  wa¬ 
ter  in  terms  of  lOPs.  Gordon  et  aL^  give  R  as 

R=fb,/{a  +  b,),  (1) 

where  6^  is  the  volume  backscattering  coefficient  and 
a  is  the  total  volume  absorption  coefficient,  which  is 
composed  of  the  sum  of  a^(X),  a^(X),  a^(X),  and  Up^CX), 
the  absorption  by  CDOM,  detritus,  pure  water,  and 
phytoplankton,  respectively;  see  Ref  1.  Factor/* is  a 
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constant  of  proportionality  whose  magnitude  de¬ 
pends  on  the  definition  of  reflectance.  For  reflec¬ 
tance  of  plane  irradiance  f  is  —0.33  and  R  is  defined 
as  EJEa~,  the  ratio  of  upwelling  irradiance  E^  to 
downwelUng  irradiance  Ea~,  where  both  quantities 
in  this  ratio  are  measured  just  beneath  the  water 
surface.  Equation  (1)  provides  the  accepted  wave¬ 
length  dependence  for  variety  of  other  expressions  for 
reflectance  from  ocean  water.  Thus  the  in-water  re¬ 
flectance  Rrsw,  defined  as  LJEf  by  Jerome  et 
is  also  given  in  terms  of  the  wavelength  dependence 
of  hj(p,  -t-  65)  according  to 

Rrsw  =  {flQ)hJ{a  +  b^),  (2) 

where  L„  is  the  upwelling  radiance  and  Q  =  EJL^^  is 
a  factor  given  by  Morel  and  Prieurs  that  accounts  for 
the  bidirectionality  of  the  Rrsw.  In  a  detailed  inves¬ 
tigation  of  the  properties  of  Q,  Morel  et  al.^  show  that 
Q  varies  with  wavelength  in  a  complex  manner. 
This  fact  would  limit  the  usefulness  of  Eq.  (2).  How¬ 
ever,  it  appears  that  the  quotient //Q  -  0.1  remains 
nearly  constant,  independently  of  wavelength  and  so¬ 
lar  angle,  for  protocol  observation  of  the  Rrsw.i 

The  most  commonly  used  expression  for  the  reflec- 
tonce  from  ocean  water  is  Rrs.  It  is  the  all- 
important  reflectance  in  optical  monitoring  of  oceans 
by  satellites.  It  is  defined  as  LJE^,  the  ratio  of 
water-leaving  radiance  and  irradiance  E^  just 
above  the  water  surface.  By  accounting  for  the 
transmittance  of  light  at  the  air— water  and  water-air 
interfaces,  denoted  tt' ,  and  taking  into  account  the 
spread  in  the  solid  angle  as  emerges  from  water  to 
become  L„,  we  obtain 

Rrs  =  {(/7Q)(tr/n/)}6fc/(a-l-6»)  [sr"'],  (3) 

where  is  the  relative  index  of  refraction  for  the 
air-water  interface  [for  Eq.  (3)  we  assume  flat  surface 
conchtions,  and  the  equation  is  not  exact;  see  Ref.  1]. 
Again,  the  coefficient  in  braces  in  Eq.  (3)  is  deemed 
nearly  independent  of  X. 

By  referencing  all  the  above  expressions  for  reflec¬ 
tance  to  plane  irradiance  we  in  essence  disregard 
the  angular  distribution  of  the  incident  light.  How¬ 
ever,  scattering  from  marine  particles  is  highly  direc¬ 
tional;  therefore  reflectance  from  ocean  water  is  also 
a  directional  quantity  whose  magnitude  and  spectral 
shape  depend  on  the  angular  distribution  of  the  in¬ 
cident  light.  Thus  referencing  Rrs(\)  to  E^  makes 
Rrs(X)  inherently  ambiguous  because  in  principle  we 
can  obtain  the  same  E^(X)  for  more  than  one  intensity 
of  incident  light.  F  or  instance,  two  different  angular 
distributions  of  daylight  illuminating  the  same  water 
mass  can  give  the  same  E^  but  different  Rrs  and  vice 
versa.  One  could  argue  that  such  differences  are 
neghgible  in  practice  for  good  protocol  observational 
conditions  that  specify  the  observational  angles  for 
Rrs  relative  to  the  Sim.  Perhaps  this  is  true,  but  we 
shall  see  that  there  is  another,  more  tangible,  diffi¬ 
culty  when  it  comes  to  the  question  of  uniqueness  of 
Rrs(X). 

Inasmuch  as  the  combination  of  factors  in  front  of 


the  term  bf,/ (a  +  b^)  in  Eqs.  (1)— (3)  is  usually  treated 
as  a  constant  of  proportionality,  we  tacitly  imply  that 
the  wavelength  dependence  of  Rrs(X)  comes  from  the 
wavelength  dependence  of  bjia  +  65).  Ladner  et 
alJ  used  in  situ  measurements  of  b^ik)  together  with 
measurements  of  a{\)  to  show  that  bjia  +  bf,)  does 
not  correlate  with  Rrs(X)  in  the  visible  region  of  the 
spectrum.  Some  of  the  discrepancy  presented  by 
Ladner  et  aU  may  come  from  the  lack  of  accuracy  in 
measurements  of  bi,{\).  Determination  of  6i(X)  is 
vulnerable  to  error  when  fej(X)  is  measured  at  one  or 
a  few  angles  that  avoid  the  anomalous  backscattering 
at  180®;  see  Ref  8. 

Sydor®  provides  an  alternative  solution  for  Rrs(X) 
that  is  not  based  on  6j(X).  Using  statistical  consid¬ 
eration  of  multiple  scattering,  Sydor®  gives  Rrs  at  any 
X  in  terms  of  the  average  number  of  scatters  per 
photon,  b/a.  The  full  relationship  between  Rrs  and 
b/a  is  given  by 

Rrs/(1  -  2iTRrs)  —  C{,6/a,  (4) 

where  Cf,  has  the  value  0.001  ±  0.0002  for  waters 
that  obey  Petzold^®  volume  scattering  function  for 
coastal  waters.  The  term  (1  —  2'7rRrs)  in  relation  (4) 
comes  from  the  correction  for  the  number  of  scatters 
per  photon  in  a  semi-infinite  medium.  For  an  infi¬ 
nite  medium  the  number  of  scatters  per  photon  is  ra  = 
(6/ a).  The  number  of  scatters  per  photon  in  a  semi¬ 
infinite  medium  reduces  to  {b/a){l  -  2-irRrs)  because 
of  reflectance.  Photons  can  escape  the  semi-infinite 
medium  before  they  scatter  an  average  of  n  times. 
For  Rrs  <Si  l/(2'n-)  the  left-hand  side  of  relation  (4) 
becomes  Rrs,  and  relation  (4)  simplifies  to  Rrs  = 
Cb{b/a).^^  Relation  (4)  provides  an  alternative  solu¬ 
tion  to  Eq.  (3)  because  it  gives  Rrs  in  terms  of  the 
dimensioifless  angle-independent  ratio  b/a  that  is 
easier  to  verify  experimentally.  Measurement  of 
a(X)  and  6(X)  is  routine  in  ocean  optics.  Relation  (4) 
differs  fundamentally  from  Eq.  (3)  because,  in  the 
former,  6(X)  and  66(X)  need  not  have  the  same  wave¬ 
length  dependence.  By  using  the  dependence  of  Rrs 
on  b/a,  Sydor  et  alM  showed  that  taking  Rrs  =» 
Cbib/a)  provides  good  estimates  of  lOPs  from  Rrs, 
provided  that  Rrs  is  analyzed  in  a  sequential  manner 
that  starts  vrith  the  determination  of  6(X)  from  the 
spectral  region  where  dominates  the  absorption. 
Sydor  et  al.^^  used  experimental  data  for  the  waters 
of  Lake  Superior  and  the  Mississippi  Sound  to  show 
that  the  proportionality  between  Rrsw  and  b/a  is 
compatible  with  the  pol3momial  expression  for  Rrsw 
presented  by  Jerome  et  al.*  In  general,  relation  (4) 
states  that  Rrs  is  not  a  linear  function  of  lOPs. 

The  differences  in  solutions  such  as  Eq.  (3)  and 
relation  (4)  are  paramount  in  discussion  of  unique¬ 
ness.  Even  if  we  assume  that  bjb  is  nearly  inde¬ 
pendent  of  wavelength,  Eq.  (3)  and  relation  (4)  differ 
mathematically.  Thus,  if  both  solutions  provide 
nearly  the  same  Rrs(X),  then  Rrs(X)  cannot  be  truly 
unique.  Sydor  et  have  already  shown  that  more 

than  one  set  of  6(X)  and  a(X)  can  produce  similar 
Rrs(X)  in  the  visible  region  of  the  spectrum,  where 
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Fig.  1.  Spectral  absorption  coefficients  derived  from  relation  (4) 
applied  in  sequential  fit  to  Rrs(X)  off  the  coast  of  South  Afnca 
before  the  onset  of  ph5iioplankton  bloom.  1,  2,  3,  Predicted  a^, 
Oph,  respectively.  Rrs(X)  was  traced  from  enlarged  figures  shown 
by  Roesler  and  Boss.2  Solid  and  filled  triangles  and  circles,  Uy  and 
Op,  respectively,  taken  from  Roesler  and  Boss.  The  wavelength 
dependence  of  the  surrogate  aph*(X)  used  in  the  sequential  fit  was 
based  on  aph(X)  measured  in  the  waters  of  the  Mississippi  Sound 
(Ref.  11).  The  magnitude  of  ap^fX)  shown  here  corresponds 
roughly  to  0.038aph*(X).  No  attempt  was  made  to  back  out  the 
actual  apii(X)  spectrum  off  the  coast  of  South  Africa.  We  assumed 
that  a^(X)  dominated  a(\)  for  X  >  650  nm.  This  assumption  is 
justified  because  at  650  nm  a^gso  ~  0*35  m“^,  which  is  far  greater 
than  the  combined  values  of  a^,  Uph,  and  at  650  nm  reported  by 
Roesler  and  Boss. 


Rrs(X.)  depends  on  some  nine  variable  parameters. 
We  test  their  result  further  by  fitting  Rrs(X)  with 
measured  lOPs,  using  both  Eq.  (3)  and  relation  (4). 

2.  Experimental  Results 

By  using  neural  networks,  Roesler  and  Boss^  em¬ 
ployed  a  version  of  Eq.  (1)  modified  in  terms  of  ex¬ 
tinction  coefficient  c(X)  to  extract  lOPs  fi*om  Rrs(X) 
and  R(X).  They  made  the  assumption  that  bjb  is  a 
weak  function  of  X  and  minimized  the  ambiguity  in 
Rrs(X)  by  procedural  methods  in  their  neural  net¬ 
works  technique.  The  use  of  neural  networks  in 
analysis  of  Rrs(X)  is  novel,  but  in  fact  it  is  still  a 
least-squares  fit.  We  performed  a  similar  task  on 
the  data  of  Roesler  and  Boss^  but  employed  a  least- 
squares  fit  to  Rrs(X),  using  relation  (4)  and  employing 
a  sequential  parameter  fit  that  starts  with  the  deter¬ 
mination  of  6(X)  in  the  spectral  region  where  a(X)  is 
dominated  by  a^^CX). 

First  we  considered  the  Rrs  data  of  Roesler  and 
Boss^  taken  off  the  coast  of  South  Africa.  Using  a 
simplified  form  of  relation  (4)  we  estimated  the  mag¬ 
nitude  of  6730  from  the  slope  of  Rrs(X)  versus  \CJ 
(Xa  )]  for  X  >  650  (nm),  where  we  took  =  0.0011 
sr“^  and  assumed  that  6(X)  behaves  as  ~(1/X).  The 
result  yielded  6(X)  ~  0.4(730/X)  Subsequently 

we  determined  the  components  of  a(X)  by  following 
the  procedure  described  by  Sydor  et  al.^^  The  re¬ 
sults  of  this  procedure  are  in  close  agreement  with 
the  experimental  data  for  Rrs(X),  as  shown  in  Fig.  1. 


Wavelength  X(nm) 

Fig.  2.  Filled  squares,  Rrs(X)  off  the  coast  of  South  Africa  taken 
from  Ref.  2.  Solid  curve,  sequential  fit  to  Rrs(X)  based  on  relation 
(4)  with  Cfe  -  0.0011  and  6(X)  =  0.39(730/X)^  \ 


The  resultant  fit  to  the  Rrs  of  Roesler  and  Boss^  is 
shown  in  Fig.  2.  The  correlation  for  the  fit  in  Fig.  2 
gave  =  0.99,  as  demonstrated  in  Fig.  3.  Thus  we 
obtained  a  result  equivalent  to  that  of  Roesler  and 
Boss2  by  using  a  different  fitting  routine,  different 
lOP  parameters,  different  assumptions,  and  different 
equations,  confirming  that  the  relationship  between 
Rrs  and  lOPs  is  ambiguous. 

To  examine  the  root  cause  of  this  ambiguity  fur¬ 
ther,  we  used  the  data  of  Roesler  and  Boss^  for  i?(X) 
taken  off  the  coast  of  Oregon.  The  use  of  i2(X)  rather 
than  Rrs(X)  is  important  because  i?(X)  is  presumably 
devoid  of  the  inaccuracies  that  are  due  to  atmo¬ 
spheric  and  surface  corrections.  By  using  the  se¬ 
quential  least-squares  fit  to  i2(X)  we  again  estimated 
6(X)  from  a  plot  ofi?(X)  versus  (C^ '/Xu^^),  assuming  in 
this  case  that  aj^(X)  dominates  a(X)  for  X  >  560  nm. 
For  reflectance  i?,  C^,'  /tV)Cij  ~  0.0065.  The 

procedure  5delded  fe{X)  ---  0.1(730/X)  m“^.  This  mag¬ 
nitude  of  6(X)  indicates  that  the  scattering  in  the 
waters  off  the  Oregon  Coast  is  dominated  by  sus¬ 
pended  particles.  Molecular  scattering  is  60  times 


Reported  Rrs  fsr’) 

Fig.  3.  Solid  line,  =  0.99  correlation  between  sequentially 
simulated  Rrs(X)  and  reported  Rrs(X)  (filled  squares  from  Fig.  2). 
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Wavelength  X(nm) 

Fig.  4.  Pilled  squares,  reflectance  R{\)  off  the  coast  of  Oregon. 
The  data  were  traced  from  Roesler  and  Boss.2  Close  fits  to  this 
J?(\)  can  be  obtained  by  use  of  real  and  false  sets  of  lOP  parame¬ 
ters.  Dotted  curve,  first-order  sequential  fit  without  any  itera¬ 
tion.  Solid  curve,  sequential  parameter  fit  by  iteration  to  find  the 
best  fit.  Dashed  curve,  example  of  an  ad  hoc  parameter  fit.  The 
contrived  ad  hoc  fit  departs  from  the  experimental  R{\)  at  longer 
wavelengths  where  dominates  the  absorption  and  thereby  lim¬ 
its  the  possibility  of  a  false  fit.  Curve  1,  fit  from  use  of 0.00656 /a^ 
alone,  indicating  that  in  this  case  1/ (Xa^y)  determines  the  spectral 
shape  of  R(\)  for  X  >  500  nm. 


loweri  than  our  6(X)  ~  0.1(730/X)  As  a  result, 
in  appl3Tng  the  sequential  fit  to  R{k)  we  ignored  mo¬ 
lecular  scattering.  On  the  other  hand,  Roesler  and 
Boss®  included  molecular  scattering  in  their  neural 
fitting  routine  for  J2(X);  thus  they  increased  the  num¬ 
ber  of  variable  parameters  even  further  than  that 
used  in  our  least-squares  sequential  fit.  Figure  4 
shows  the  least-squares  fit  to  R{k)  off  the  Oregon 
Coast.  The  correlation  for  J?(X)  was  ~  0.99,  com¬ 
parable  with  the  result  of  Roesler  and  Boss®  but  with¬ 
out  the  use  of  the  extra  parameters  that  are  due  to 
molecular  scattering. 

The  results  in  Fig.  4  demonstrate  the  ambiguity  of 
reflectance  fi'om  ocean  water  in  the  visible  region  of 
the  spectnun  and  simply  point  to  the  fact  that  any 
broad  fimction  of  X  such  as  Rrs(X)  or  R(k)  can  be 
approximated  closely  by  more  than  one  combination 
of  broad  functions  such  as  6(X),  ay(X),  and  a^(X)  that 
together  have  several  adjustable  parameters;  the 
greater  the  number  of  parameters,  the  greater  the 
ambiguity.  To  demonstrate  this  ambiguity  caused 
by  overparameterization  of  i?(X)  we  produced  an  ad 
fioc  fit  to  R(X)  by  neglecting  aph(X)  altogether  and 
varying  at  will  the  magnitude  and  the  wavelength 
dependence  of  b{k)  and  a(k)  without  Cpi,  until  we 
came  up  with  a  close  fit  to  the  experimental  R{k). 
The  ad  hoc  fit  is  also  shown  in  Fig.  4.  We  can  see 
that  it  approximates  the  measured  R{k)  with  ~ 
0.98  or  better.  The  scattering  and  the  absorption 
coefficients  for  the  sequential  and  the  ad  hoc  fits  are 
shown  in  Figs.  5  and  6,  respectively.  The  ad  hoc  b(k) 
in  Fig.  5  is  deemed  unlikely  because  its  magnitude 
increases  with  the  wavelength,  whereas  b(k)  usually 


Wavelength  X  (nm) 

Fig.  5.  Curve  1,  6(X)  determined  from  R{\)  versus  C^fa^,  in  the 
X  >  560  nm  region  of  jR(X)  shown  by  the  filled  squares  in  Fig.  4. 
Curve  2,  6(X)  for  the  false  ad  hoc  fit  to  R{\).  6(X)  obtained  in  the 
ad  hoc  fit  is  deemed  unlikely,  as  its  magnitude  increases  with  X. 
Nonetheless,  the  ad  hoc  6(X)  produces  a  close  fit  to  R(\)  as  long  as 
a^{\)  and  a^(X)  are  adjusted  accordingly. 

decreases  as  ~(1/X).  Yet  the  unlikely  b{k)  provides 
an  excellent  fit  to  R(k)  for  X  <  500  nm  as  long  as  we 
adjust  a(X)  accordingly.  Basically,  any  unrestricted 
fit  to  Rrs(X)  would  lead  to  a  similar  result,  because 
any  given  Rrs(X)  does  not  translate  to  a  unique  com¬ 
bination  of  lOPs  in  ocean  water;  see  Ref  1.  Clearly 
the  situation  is  not  hopeless  because  we  can  restrict 
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Fig.  6.  The  solid  trace  and  curve  2  show  a^(X)  from  the  sequential 
and  the  ad  hoc  fits,  respectively,  toi?(X)  off  the  Oregon  Coast  (filled 
squares  in  Fig.  4).  Curve  3,  a^(X),  the  detrital  component  of  a^CX). 
Most  of  the  absorption  by  suspended  particles  off  the  Coast  of 
Oregon  is  attributable  to  OphCX),  whose  magnitude  here  was 
0.0035aph*(X)  of  the  waters  of  the  Mississippi  Sound.  The  lone 
upward-pointing  triangle  shows  an  experimental  point  for  Op  re¬ 
ported  by  Roesler  and  Boss.^  Curves  A  and  B,  a^(X)  for  the  se¬ 
quential  and  the  ad  hoc  fits,  respectively.  The  downward¬ 
pointing  triangle  shows  the  lone  a^  point  reported  by  Roesler  and 
Boss. 2  We  assumed  that  dominates  the  absorption  in  the  tail 
ofR(X)  for  X  >  560  nm.  Curve  1  of  Fig.  4  shows  that  this  assump¬ 
tion  was  justified.  The  pure-water  tail,  where  a^(K)  —  a(X),  shifts 
toward  the  shorter  wavelengths  when  the  concentration  of  sus¬ 
pended  particles  and  CDOM  is  low. 
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and  prioritize  the  search  for  reasonable  lOPs  by 
knowing  the  general  properties  of  lOPs  for  any  area 
ahead  of  the  time.  Nonetheless,  the  fact  that  any 
given  Rrs(X)  does  not  represent  a  unique  set  of  lOPs 
places  a  limit  on  the  accuracy  of  the  remote-sensing 
determination  of  lOPs,  especially  aph(X),  whose  a  pri¬ 
ori  wavelength  dependence  is  always  uncertain. 

The  upshot  of  the  above  presentation  is  that  we 
need  some  prior  knowledge  of  the  wavelength  depen¬ 
dence  of  lOPs  when  we  try  to  determine  lOPs  re¬ 
motely.  For  instance,  any  fitting  routine  that  takes 
advantage  of  the  fact  that  l/[Xnj„(X.)]  determines  the 
spectral  shape  of  Rrs(X)  in  a  long-wavelength  tail  is 
likely  to  select  a  realistic  value  for  6(\)  and  subse¬ 
quently  determine  a  realistic  a(X)  relative  to  the 
known  magnitude  ofajj,(X).^^  Consider,  for  example, 
R(X)  of  off  the  Coast  of  Oregon.  Half  of  its  spectral 
shape,  the  region  X  >  500  nm,  is  determined  solely  by 
l/(Xa^),  as  demonstrated  by  curve  1  in  Fig.  4.  In 
general,  we  can  provide  reasonable  approximations 
to  Rrs(X)  or  R(k)  based  on  the  expected  exponential 
dependence  of  c^(X)  and  Oaik)  and  the  ~1/X  depen¬ 
dence  of  6(X).  Subsequently  we  can  back  out  api,(X) 
that  provides  a  perfect  fit  to  Rrs(X)  as  shown  in  Ref. 
11.  However,  we  can  never  be  sure  of  the  accuracy 
of  OpiiCX)  because  we  are  never  sure  of  the  svrface 
corrections  for  Rrs(X)  and  the  inherent  ambiguity  in 
Rrs  that  is  due  to  the  uncertainty  in  the  angular 
distribution  of  the  incident  light. 

Much  of  the  ambiguity  that  is  attributable  to  over¬ 
parameterization  comes  from  the  fact  that  and 
have  similar  overlapping  exponential  dependence  on 
X.^i  Thus  the  relative  contributions  from  Uy  and 
to  the  total  a(X)  can  produce  a  variety  of  close  approx¬ 
imations  to  any  given  Rrs  spectrum.  Clearly,  no 
matter  what  fitting  routines  we  use,  the  ambiguity 
produced  by  the  similarity  of  the  wavelength  depen¬ 
dence  of  OyfX)  and  a^ik)  could  never  be  resolved  if 
Uyik)  and  a^ik)  had  exactly  the  SEune  wavelength  de¬ 
pendence.  The  overparameterization  vanishes  in 
the  tail  of  Rrs(X)  because  Oy  and  drop  off  exponen¬ 
tially  while  increases  rapidly  as  X  ^  700  nm.  For 

X  >  750  nm,  all  absorption  coefficients  become  neg¬ 
ligible  except  for  a^;  thus  we  can  use  this  region  of  the 
spectrum  to  isolate  the  dependence  of  Rrs(X)  on  6(X). 

To  examine  the  discrepancy  in  Rrs(X)  represented 
by  Eq.  (3)  versus  Rrs  expressed  in  terms  of  b/ a,  we 
examined  Rrs  for  Pearl  River,  Miss.,  and  the  Gulf  of 
Mexico.  Figure  7  shows  Rrs(X)  for  the  turbid  waters 
of  the  Pearl  River.  It  also  shows  the  least-squares  fit 
fi*om  relation  (4)  and  the  expected  Rrs(X)  using  in  situ 
measimements  of  6<,(X),  6(X),  aJk),  aa(k),  and  Oph(X)  in 
Eq.  (3)  and  relation  (4).  All  magnitudes  of  Rrs  in 
Fig.  7  were  normahzed  to  the  measured  Rrs  at  620 
nm  to  permit  comparison  of  the  wavelength  depen¬ 
dence  predicted  by  each  equation.  The  measured 
Emd  simulated  lOPs  of  Pearl  River  waters  are  shown 
in  Figs.  8  and  9.  If  we  take  ^^(X)  measured  by  use  of 
a  Hydroscat  instrument  at  face  value,  the  result 
shown  in  Fig.  7  demonstrates  that  Eq.  (3)  does  not 
produce  a  close  fit  to  Rrs(X).  The  shortcoming  of  Eq. 
(3)  applied  to  measuring  the  highly  turbid  Pearl 
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Fig,  7.  Wavelength  dependence  of  Rrs(X)  for  turbid  waters  of 
Pearl  River,  Miss.  Solid  and  dashed  curves,  two  (independently) 
surface  corrected  values  of  Rrs(X)  for  Pearl  River.  The  triangles 
show  that  ^^/(a  +  h^)  does  not  follow  the  measured  Rrs(X);  the 
filled  squares  show  that  (6/a)(l  —  2iTRrs)  follows  the  spectral 
shape  of  Rrs(X)  quite  closely.  To  compare  their  spectral  depen¬ 
dence  we  normalized  both  cases  to  the  measured  Rrs  at  620  nm. 
The  dotted  curve  shows  that  the  sequential  fit  to  Rrs(X)  has  the 
correct  magnitude  and  displays  the  fine  spectral  features  exhibited 
in  the  measured  Rrs(X). 


River  is  not  surprising.  Whitlock  et  suggested 
that  Eq.  (3)  fails  for  turbid  waters.  Ladner  et  alJ 
demonstrated  a  similar  result  in  their  Fig.  5.  In 
fact,  their  data  show'^  that  Rrs  does  not  correlate  well 
with  bjia  +  65)  for  X  <  650  nm  but  that  the  corre¬ 
lation  improves  markedly  for  X  >  676  nm,  where  the 
wavelength  dependence  of  fc^CX)  and  6(X)  becomes 
similar,  as  indicated  in  Fig.  8.  The  similarity  of  the 
wavelength  dependence  of  6(X)  and  fe^CX)  for  X  >  676 
nm  could  be  attributed  to  the  fact  that  may  affect 
6^  (owing  to  the  electromagnetic  boundary  conditions 
on  the  illuminated  side  of  the  particles)  more  than  it 


0 1 — ^ — I — ^ — I — ^ — I — . — I — ^ ^ ^ — > — ■ — ' 
400  450  500  550  600  650  700  750  800 


Wavelength  X(nm) 

Fig.  8.  Filled  squares,  b (X)  for  Pearl  River  determined  with  an  ac9 
instrument.  Stars,  6e,(X)  measured  with  Hydroscat  (multiplied 
here  by  a  factor  of  14  to  permit  its  wavelength  dependence  to  be 
compared  with  that  of  6(X)).  Solid  curve,  6(X)  obtained  from  a 
sequential  fit  by  use  of  relation  (4)  and  =  0.0012  sr 
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Fig.  9.  Downward-pointing  triangles  a(\)  (without  aj  for  the 
waters  of  Pearl  River,  Miss.  Upward-pointing  triangles,  the  sum 
of  a^,(X)  measured  with  a  0.045-pm  pore  filter  at  the  intake  of  the 
ac9  meter  plus  a^CX)  determined  fii'om  a  filter-pad  transmission 
measurement.  Solid  curve,  a(X)  obtained  fi-om  the  sequential  fit 
fi-om  relation  (4). 


affects  the  total  6.  Thus,  when  a^(\)  ->•  0  as  \  700 

nm,  bi,  and  b  attain  a  similar  wavelength  dependence 
that  depends  only  on  the  particle  size  distribution 
and  on  the  average  index  of  refraction;  then  the  cor¬ 
relation  between  Rrs  and  bb/ia  6^,)  improves. 

By  comparing  the  values  of  Rrs  predicted  by  Eq.  (3) 
and  relation  (4)  we  can  see  from  Fig.  7  that  relation 
(4)  produces  a  closer  fit  to  the  measured  Rrs(\)  be¬ 
cause  it  accoimts  for  all  scattering.  What  can  we 
conclude  from  such  results?  The  imderlying  laws  of 
physics  are  not  ambiguous  and  should  produce  only 
one  correct  result.  If  diffuse  reflectance  were  refer¬ 
enced  to  the  incident  radiance  rather  tham  irradiance, 
Rrs  would  be  defined  uniquely,  but  the  problem  of 
overspecification  would  remain.  Terms  such  as  f/Q 
^e  designed  to  accovmt  for  the  distribution  of  the 
incident  radiemce  and  multiple  scattering.  How- 
eyer,  f/Q  introduces  its  own  ambiguity  because  addi¬ 
tional  variables  tend  to  disguise  the  inadequacy  of 
Eq.  (3)  that  arises  from  the  fact  that  Rrs  cannot  be 
described  solely  in  terms  of  the  plane-wave  single 
backscattering  process  that  is  implied  by  use  of  6^  in 
Eq.  (3).  For  extended  illumination  and  multiple 
scattering,  forward  scattering  also  contributes  to  the 
spectral  dependence  of  Rrs(\).  Indeed,  photons  ob¬ 
served  in  Rrs  come  from  variety  of  combinations  of 
forward  and  backward  scattering  whose  combined 
average  wavelength  dependence  governs  the  wave¬ 
length  dependence  of  Rrs(X).  In  multiple  scattering 
each  reflected  photon  samples  a  variety  of  particles  in 
a  process  that  is  fundamentally  different  from  that 
represented  by  b^,,  where  by  definition  each  photon 
samples  one  particle  only.  Thus  in  multiple  scatter¬ 
ing  the  sampling  of  the  optical  properties  of  sus¬ 
pended  particles  is  interrelated  from  one  scatter  to 
the  next.  As  a  result,  Rrs  from  multiple  scattering  is 
not  additive,  as  implied  by  single-scattering  equation 
(3)  for  which  each  photon  scatters  only  once  and  sam- 


Fig.  10.  In  the  linear  limit  of  Rrs  1  and  constant  bjb  <  0.03, 
Eq.  (3)  and  relation  (4)  give  comparable  results  for  Rrs(X.). 
Downward-pointing  triangles,  Rrs(\)  obtained  from  Eq.  (3)  by  Hy¬ 
droscat  6(,(X)  and  ac9  measurement  of  a(X).  Upward-pointing  tri¬ 
angles,  Rrs(X)  predicted  by  relation  (4)  from  ac9-determined  6(X) 
and  a(X).  Solid  curve,  measured  Rrs(X). 


pies  the  optical  properties  one  particle  at  a  time  in¬ 
dependently  of  other  particles.  Rrs  from  multiple 
scattering  is  not  a  linear  function  of  lOPs.  Relation 
(4)  includes  multiple  scattering  and  photon  turn¬ 
around  through  computation  of  the  magnitude  of  C*.9 
It  states  that  average  Rrs(X,)  is  proportional  to  (6/ 
®)(1  ~  2'ji’Rrs),  the  average  munber  of  scatters  per 
photon  in  a  semi-infinite  medium.  Thus  relation  (4) 
states  that  Rrs  =  0  only  if  6  =  0.  Equation  (3)  does 
not  account  for  photon  turnaround  and  states  that 
Rrs  ^  0  if  6ft  ->  0  even  if  6  is  finite.  This  is  clearly 
erroneous.  Equation  (3)  and  relation  (4)  approach 
the  same  result  in  the  linear  limit  of  Ifrs  <sc  1  and 
constant  6ft/6  <  0.03,  as  demonstrated  experimen¬ 
tally  by  Rrs(X)  in  Fig.  10.  In  the  linear  limit,  bino¬ 
mial  expansion  of  Eq.  (3)  and  relation  (4)  shows  that 
they  differ  largely  in  the  magnitude  of  the  propor¬ 
tionality  constant. 

3.  Need  for  Hyperspectral  Data 

We  can  resolve  the  ambiguity  from  overparameter¬ 
ization  by  identifying  the  region  of  the  spectrum 
where  dominates  the  absorption.  However,  the 
pure-water  region  shifts  toward  the  shorter  wave¬ 
lengths,  depending  on  the  concentration  of  suspended 
particles  and  CDOM.  To  monitor  a  wide  variety  of 
oce^  waters  one  needs  Rrs(\)  data  that  cover  the 
entire  380-950-nm  spectrum  with  a  sufficient  num¬ 
ber  of  bands  that  the  location  of  the  shifting  spectral 
region  where  aj\)  ~  a(X)  can  be  identified.  The 
upper  wavelength  limit  at  950  nm  is  needed  as  the 
point  of  reference  where  Rrs  becomes  negligible  be¬ 
cause  Cu,  30  m“^;  i.e.,  where  we  can  set  the  mag¬ 
nitude  of  the  atmospherically  corrected  Rrs  at  zero 
and  use  it  to  estimate  the  magnitude  of  the  rough- 
surface  reflectance.il  Similarly,  Shybanovi^  points 
out  that  the  region  near  380  nm  is  useful  for  estimat¬ 
ing  the  magnitude  of  Oy. 
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4.  Conclusions 

Rrs(X)  does  not  correspond  to  a  unique  set  of  lOPs. 
Values  of  lOP  close  to  those  of  the  actual  lOPs  can  be 
obtained  from  spectral  analysis  of  Rrs(X),  provided 
that  we  impose  realistic  limits  on  the  wavelength 
dependence  of  6(X),  a^CX),  and  a^(X).  Usually  fe(X) 
has  a  smooth  monotonic  ^1/X  dependence,  imlike 
6^(X).  Thus,  by  using  relation  (4)  under  the  condi¬ 
tion  that  Rrs  «  1  and  a  -  a^,  we  can  estimate  the 
magnitude  of  6(X)  from  plots  of  Rrs(X)  versus  C^/ 
(ka^).  Cl,  varies  by  20%  with  particle  concentration 
but  appears  constant  as  a  function  of  X,  and  its  aver¬ 
age  magnitude,  0.001  sr“^,  appears  to  hold  for  a  va¬ 
riety  of  coastal  waters.  However,  it  is  unlikely  that 
Cl,  is  global  outside  the  hmits  0.01  <  bjh  <  0.04  set 
by  Jerome  et  al^  It  is  also  unlikely  that  Eq.  (3)  is 
global  for  open  ocean  waters  where  CDOM  is  present. 

We  have  made  use  of  data  presented  by  Roesler 
and  Boss^  and  by  Ladner  et  al!^  published  in  the 
proceedings  of  the  Ocean  Optics  XVI  Conference  in 
Santa  Fe,  N.Mex.  on  18-23  November  2002.  We  are 
indebted  to  those  investigators  for  sharing  their  re¬ 
sults.  We  are  also  grateful  to  our  colleague  A1  Wei- 
dermann  for  personal  commxmications  and  for 
providing  us  with  data  on  the  spectral  dependence  of 
the  scattering  phase  function  measured  by  Misha  E. 
Lee  and  Eugene  B.  Shybanov. 
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